Generation and Structure of Solitary Rossby Vortices in Rotating Fluids 
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The formation of zonal flows and vortices in the generalized Charney-Hasegawa-Mima equation is 
studied. We focus on the regime when the size of structures is comparable to or larger than the 
deformation (Rossby) radius. Numerical simulations show the formation of anticyclonic vortices in 
unstable shear flows and ring-like vortices with quiescent cores and vorticity concentrated in a ring. 
Physical mechanisms that lead to these phenomena and their relevance to turbulence in planetary 
atmospheres are discussed. 



PACS numbers: 47.27.-i, 92.60 Ek 

Fluid motion subject to strong rotation, stratification, 
or action of a magnetic field can often be considered as 
quasi-two-dimensional. One of the major features of 2D 
flows is their tendency to self-organization which reveals 
itself by spontaneous generation of coherent structures 
(e.g. vortices, jets) that dominate the large-scale mo- 
tion. The dynamics of these structures becomes more 
complicated when the inverse energy cascade interferes 
with the characteristic spatial scale determined by the 
external field. The presence of the Coriolis force due to 
rotation (/3-effect) leads to a limitation on vortex sizes in 
the meridional direction by the Rhines length [jlj , while in 
the strong turbulence regime when waves are neglected, 
the influence of the finite deformation (Rossby) radius re- 
sults in the emergence of "shielded" vortices, which can 
form a quasicrystalline structure . 

In this Communication we extend previous studies to 
include both wave and finite deformation radius effects, 
so that the model becomes more relevant to geophysics 
and plasma physics applications and reveals new physical 
effects. It is known from observations that, in planetary 
atmospheres, there exist global-scale circulations in the 
form of longitudinal zonal flows with embedded long-lived 
vortices. Indeed, this effect is significant in the atmo- 
spheres of giant planets and in the Earth's oceans. The 
mechanism of the formation of zonal flows on the surface 
of giant planets is a separate question which is not well 
understood. Two distinct approaches to this problem 
are based respectively on 3D thermal convection H and 
on the inverse energy cascade in 2D turbulence. Very 
recent data obtained by the Galileo probe during the 
unique first-ever in situ measurements in the atmosphere 
of Jupiter g seem to be more in favor of the 3D convec- 
tion concept: the observations showed stronger than pre- 
viously assumed winds (up to 200 m/s) and turbulence 
in the upper layer of the Jovian atmosphere. This indi- 
cates that the origin of Jupiter's winds and circulation 
patterns is probably heat escaping from Jupiter's deep 
interior. However, the eddies we are interested in here 
are believed to be quasi-2D structures originated from 
the instability of zonal flows and confined to a shallow 
atmospheric layer. 

There have been several studies devoted to the com- 
bined impact of finite deformation radius and the /3-effect 



|3|,|6|. However, it is difficult to detect any particular new 
physical effects, since the scales overlap thus making the 
whole picture vague. If the rotation is strong enough, 
zonal flows destroy coherent vortices. In the opposite 
case, when the deformation radius is smaller than the 
Rhines length, the formation of "shielded" weakly inter- 
acting vortices significantly reduces the inverse cascade, 
so that zonal flows do not form. We will show how this 
difficulty can be partly overcome if the simplistic model 
is modified to take into account additional effects that 
take place when scales of the order or larger than the 
deformation radius are considered. 

We consider the generalized geostrophic (or Charney- 
Hasegawa-Mima) equation for Rossby- wave turbulence in 
the /3-plane approximation: 

d t {V 2 h - X 2 h) + [h, V 2 h] + [3d x h(l + h) = D + F, (1) 

where A is the ratio of a spatial scale used for nor- 
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malization to the Rossby radius Lr — (gHg) //o, g 
is the gravitational acceleration, h = (H — Hq)/Hq is 
the perturbation of the atmosphere of average depth H , 
f = fo+Py is the Coriolis parameter, [a, b] — a x b y — a y b x , 
D and F is the dissipation and forcing respectively, and 
the following dimensionlcss variables are used: x,y — > 
XLrx, XL R y, t^X 2 t/f Q , (3 -> f P/X 3 L R . Eq.(l) 
is isomorphic to the equation for drift waves in a mag- 
netized plasma where the inhomogeneity of the Coriolis 
parameter is replaced by the gradient of the electron den- 
sity (see, e.g., for details). Consequently all results can 
be transferred to the plasma case. The equation contain- 
ing the quadratic /3-term was derived for the first time in 
H and has been studied since then for various applica- 
tions The presence of the scalar nonlinearity (3hd x h 
is known to introduce cyclone/ anticyclone asymmetry in 
Eq.(l) (in a cyclone the direction of rotation coincides 
with that of the system). This term appears due to the 
perturbation of the fluid depth under the influence of the 
Rossby wave, so that the full depth H = H$ + h should 
be retained rather than the mean depth Hq. 

Studies of shear flows formed due to the strong rotation 
in /3-plane 2D models indicate that these flows are sta- 
ble as described by the Rayleigh-Kuo instability criterion 
jH)|,[ll]]. This result obviously contradicts the observed 
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coexistence of coherent long-lived vortices and zonal jets 
in the Jovian atmosphere. Following || we show that 
this contradiction can be resolved when the finite defor- 
mation radius is taken into account. The Rayleigh-Kuo 
criterion for shear flow instability in the approximation 
of a constant depth of the atmosphere (Hq — const) is 
given by d 2 U/dy 2 — (3q — 0, where U is the average hor- 
izontal flow velocity, y is the meridional coordinate, and 
(3q = df/dy. However, if there is a free surface, the 
average depth of the atmosphere Hq can change in the 
meridional direction: dH /dy / 0. This gradient is bal- 
anced by the Coriolis force gdHg/dy = —fU, and (3 is 
given by 
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Consequently, the modified instability criterion is 



d 2 U 
dy 2 



or in our dimensionless units 



d 2 U 
dy 2 



f3 -X 2 U = 



(3) 



(4) 



It is now clear that zonal flows which otherwise would be 
stable can become unstable when the finite deformation 
radius is taken into account. At the same time, the scalar 
nonlinearity in Eq.(l), as in the case of the Korteweg-de 
Vries equation, can lead to the gradual steepening of the 
initial perturbation in the longitudinal direction which 
is compensated by negative dispersion only for anticy- 
clones. One can thus expect that vortices emerging due 
to the instability of zonal flows will be mostly anticy- 
clones. 

We used the following setup for the numerical exper- 
iment to check this prediction. In 2D turbulence it has 
been convincingly demonstrated that the case in which 
long-term dynamics does not depend on initial condi- 
tions is best achieved with small-scale random forcing. 
In this case zonal flows are known to be a robust fea- 
ture of the flow evolution on the /3-plane. We solve 
Eq.(l) numerically using a pseudospectral method in a 
square domain 2n x 2tt with doubly periodic boundary 
conditions and resolution 512 x 512. To confine the 
dissipation to the smallest scales we use hyperviscosity 
D = {-l)P +1 v p V 2p {V 2 h), p = 8. We start with zero 
initial conditions, random forcing at 100 < kf < 105, and 
/3 = 100, A = 0. One can expect the formation of zonal 
flows to begin when the inverse energy cascade reaches 
the Rhines scale Lp ~ k^ 1 , kp = (f3 2 /e) 1 ^ 5 , where e 
is the energy injection rate, and in the considered case 
kp ~ 60. The flow is dominated by well-developed zonal 
flows by t — 15to, where tq ~ ir/v rms is the large eddy 
turnover time. At this moment the force is turned off 
to let the viscosity smooth out the irregularities in the 
smallest scales. At time t = 17tq one can observe quite 



regular shear flows. The horizontal velocity has a pro- 
file resembling the one on Jupiter with smooth westward 
and peaked eastward flows (Fig. la). These flows are 
stable, since U yy — /3q < for all y (Fig. lb). After 
that we "turn on" the deformation radius (or equiva- 
lently remove the "rigid lid" and allow the fluid to have 
a free surface) by taking A = 10. The choice of A is 
made such that the parameters would be close to the 
ones for Jupiter's Great Red Spot (JGRS). For Jupiter 
R/L R ~ 7 • 10 4 fcm/6 • 10 3 fcm - 10, f GRS ~ l.^lO -4 *- 1 , 
and since j3 ~ f/R, this gives (3 ~ 100 in our dimension- 
less units. The value U yy — (3q — X 2 U now changes sign for 
some y, and the westward anticyclonic (U < 0) zonal flow 
become unstable. This almost immediately (At ~ to) 
leads to the formation of vortices. They presumably ap- 
pear as cyclone-anticyclone pairs with cyclones decay- 
ing faster, so that only anticyclones survive (Fig. 2). To 
illustrate the asymmetry of the solution, in Fig. 3 we 

plot the vorticity skewness S u (t) = (uj 3 ) J (co 2 ) 3 ^ 2 , where 
ijj = V 2 /i, ( ) denotes an area average (u < corre- 
sponding to anticyclones). 

Another important feature to verify is the form of the 
vortices. As was discussed in |^Jl^| the scalar nonlin- 
earity in Eq.(l) allows soliton-like solutions. These solu- 
tions combine characteristic features of both waves and 
real vortices, in the sense that they have trapped fluid, 
so they are called "solitary vortices" . An analytically ob- 
tained solution H under the assumption h <C H predicts 
the structure of such a soliton to be given by 



h/ho = cosh 3 / 4 ( — ), a ~ L R h 



-1/2 



(5) 



Consequently, the vorticity profile w = V 2 h = ^■§^(r^p) 
is smooth, gradually decreasing (increasing) from the 
center of a vortex. Observations show, however, that, 
for example, the JGRS and intrathermoclinic vortices 
(lenses) in the Earth's ocean do not rotate as solid bod- 
ies, and their vorticity has a circumferential profile with 
a rotating perimeter and relatively quiet core. Also in 
experiments in rotating tanks |ll| with h ~ O(l), vor- 
tices appear to be more localized than predicted by the 
theoretical solution (5). Sutyrin noted jlij that much 
better agreement with experiments can be obtained if 
one uses the assumption that there is a uniform profile of 
the potential vorticity £ = V 2 /i — X 2 h = const inside the 
vortex, although the vorticity profile was never analyzed. 
Ironically, Marcus [l5| used the same assumption to crit- 
icize the soliton-like model, arguing that it is unable to 
produce ring vortices. 

We show here that the model Eq.(l) does produce 
ring vortices and we present the evidence that it is the 
scalar nonlinearity that plays a decisive role. It is impor- 
tant to emphasize that these vortices spontaneously arise 
from random initial conditions; they are not in any sense 
"trial" vortices introduced into the flow. 

In our numerical experiment without forcing we start 
with the following initial perturbation 
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Cexp(i<f>k) if 4 < k < 5 
otherwise, 

with random phases 0k € [0, 2tt). C is chosen such that 
total initial energy, E = \ J2k(k 2 + ^ 2 )l^k| 2 , is equal to 
0.5, and A = L^ 1 = 10, (3 = 20. The characteristic 
size of the initial eddies is L ~ (6 -f- 7)Lr. During the 
evolution of the system for over 50ro we observe the for- 
mation of elongated vorticity sheets, which deform into 
loops. These loops could close into ring vortices, all of 
them being anticyclones (Fig. 4). If the scalar nonlincar- 
ity is dropped from the equation, closed loops and ring 
vortices (in this case they could equally be cyclones and 
anticylones) do not form (Fig. 5a). In the case (3 = 
the cubic nonlinearity J(h,hV 2 h + i(Vft.) 2 ) should be 
retained in Eq.(l) j5|. The cyclone/ anticyclone asymme- 
try is then preserved, but no signs of ring vortices are 
found (Fig.5b). 

We do not have a theoretical explanation of the for- 
mation of ring vortices, but necessary physical conditions 
can be identified. First of all, because of the finite defor- 
mation radius Lr vortices are shielded, i.e. their inter- 
action decreases exponentially. This means that a ring 
vortex should be of the size of several Lr in order to be 
stable. Simultaneously, the Rhines length Lp and the 
Rossby radius Lr should be comparable; in the opposite 
case either zonal flows or "frozen", not mobile, vortices 
dominate the flow. Mutual action of negative dispersion 
and KdV-like scalar nonlinearity can lead to the forma- 
tion of loops. Connection of the opposite sides of a loop 
leads to the formation of a ring vortex. Thus, both fi- 
nite deformation radius and the scalar nonlinearity are 
crucial for the emergence and stability of ring vortices. 

In conclusion, we have studied the question determin- 
ing the combined influence of the /3-effect (both Rossby 
waves and KdV nonlinearity) and deformation radius on 
coherent structures. It is shown that for a certain set of 
parameters (close to those of giant planets) large-scale 
zonal flows created by small-scale forcing tend to pref- 
erentially form anticyclonic vortices. Their emergence 
can be explained by the modified Rayleigh-Kuo instabil- 
ity criterion, which takes into account the deformation 
of free surface. We have provided numerical evidence 
that for scales larger than the Rossby radius, the scalar 
nonlinearity is responsible for the formation of ring an- 
ticyclonic vortices, a feature known from observations. 
We have demonstrated how gradual complication of the 
simple model allows one to clearly reveal new physical 
effects. The complete and accurate study of these effects 
will need simulations of more complicated equations on 
the sphere (in the spirit of |l6| , |l7| ] ) . 
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FIGURE CAPTIONS 

Figure 1. Averaged zonal velocity U (a) and U yy 
as functions of latitude y. 
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Figure 2. Finite deformation radius makes zonal flows 
unstable with respect to generation of anticyclones. The 
vorticity field before (a) and after (b) introducing a finite 
deformation radius L r (whose scale is shown in the right 
picture). Darker colors correspond to the anticyclons. 

Figure 3. Vorticity skewness S u (t), X — before t = 15r 
and A = 10 afterwards. 

Figure 4- Instantaneous vorticity field at t — 40to show- 
ing the ring anticyclonic vortices, S u = —2.8 (a); cross 
section of a typical vortex (b). 

Figure 5. Vorticity fields at the same time as in Fig. 4 
when (a) the scalar nonlinearity is dropped, S u = 0; 
(b) ,9 = 0, but the term J(h,hV 2 h + \{Vh) 2 ) is added, 
5 W = -0.3. 
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Tignre 2: Finite deformai ion radius makes sonal flows unstable with respect to generation of 
anticyclones. The vorticity field before (a) and after (b) introducing a finite deformation radius 
Lr (whose scale is shown in the right picture). Darker colors correspond to anticyclones. 




FIG. 3. Vorticity skewness S u (t), A = before t = 15tq and A = 10 afterwards. 
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"Figure 4: Instantaneous vorticity field a,t i = 4Qr o showing the ring ant i cyclonic vortices. 
So, = —2.8 (a,); cross section of a. typical vortex (b). 




"Figure 5: Vcrticity fields at the same time as in Fig.4 when (a] the scalar noriliriearity is 
dropped, S„ = 0; (b) ft = 0, but the term J[k, AV 5 /t + j(VA) 3 ) is added, S„ = -0.3. 



